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ABSTRACT 

The nuclei of Seyfert 1 galaxies exhibit a range of optical polarization characteristics that can be 
understood in terms of two scattering regions producing orthogonal polarizations: an extended polar 
scattering region (PSR) and a compact equatorial scattering region (ESR), located within the circum- 
nuclear torus. Here we present NICMOS 2.0/xm imaging polarimetry of 6 "polar scattered" Seyfert 1 
(SI) galaxies, in which the PSR dominates the optical polarization. The unresolved nucleus (< 0'.'58) 
is significantly polarized in only three objects, but 5 of the 6 exhibit polarization in a 0'.'58-l'.'5 circum- 
nuclear annulus. In FairallSl and ESO323-G077, the polarization position angle at 2fim (02/im) is 
consistent with the average for the optical spectrum (^t,), implying that the nuclear polarization is 
dominated by polar scattering at both wavelengths. The same is probably true for NGC 3227. In 
both NGC 4593 and Mrk766, there is a large difference between 6'2/im and 9^ off nucleus, where polar 
scattering is expected to dominate. This may be due to contamination by interstellar polarization in 
NGC 4593, but there is no clear explanation in the case of the strongly polarized Mrk766. Lastly, in 
Mrkl239, a large change (« 60°) in 02/im between the nucleus and the annulus indicates that the 
unresolved nucleus and its immediate surroundings have different polarization states at 2/im, which we 
attribute to the ESR and PSR, respectively. A further implication is that the source of the scattered 
2/xm emission in the unresolved nucleus is the accretion disk, rather than torus hot dust emission. 

Subject headings: galaxies: nuclei - galaxies: Seyfert - infrared: galaxies - polarization - scattering - 
techniques: polarimetric 



1. INTRODUCTION 

In the Unified Model for S eyfert galaxies (jAntonuccil 
[T993I: lUrrv fc Padovanil[T995h . Seyfert Type 1 (SI) and 
Type 2 (S2) nuclei are intrinsically the same type of ob- 
ject viewed at different orientations. In the S2's, our 
direct line-of-sight to the nuclear continuum source and 
broad-line region (BLR) is blocked by a clumpy toroidal 
region of dusty molecul ar gas clouds, 011 scales of possi- 
bly j ust a few parsecs (jJaflte et al.l 120041 : iTristram et al.l 
12007*). Spectropolarimetry played a pivotal role in es- 
tablishing this picture through the detection of polar- 
ized broad-lines. These features, attributed to scat- 
tering of b road-line emission above the poles of the 
torus (e.g., lAntonucci fc Milleil 119851 ). reveal the pres- 
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ence of an otherwise obscured BLR in many S2's. The 
search for polarized broad-lines has motivate d many sub- 
sequent spectropolarimetric studies of S2's ()Tran et al.l 
[T99I I Young et al.l [19961 1^312001 . The optical polar- 
ization position angle (0„) in S2's is usually oriented per- 
pendicular to the projected radio sourc e axis and hence 
t he axis of the obscuring toroidal region (|Antonuccilll983l: 
Brin dle et al.l 11990). This is consistent with the simple 
polar scattering envisaged in the unified model, because 
scattered light is polarized perpendicular to the scatter- 
ing plane that contains the incident rays. 

In contrast, 9^ is parallel to the radio source axis in 
the majority of Si's (| Antonuccill 19831: iSmith et al.| [2002'). 
This implies scattering in a plane perpendicular to the 
system principal axis, which in turn indicates a second 
scattering region that is present in Si's but not observed 
in S2's. Spectropolarimetric studies have shown that 
Si's often exhibit distinctive structure in both 9^ and the 
percentage of optical polarization (p„) across the broad 
Ha emission- line profile ( Goodrich & Miller 1994; MartH 
IT9961 IT998t lYoung et al.l1l999t ISmith et al.ll2002D . Such 
features are naturally produced if the line emission orig- 
inates in a rotating disk (presumably the outer regions 
of the accretion disk itself) and is scattered in a com- 
pact region that is co-planar with the disk, and closely 
surrounds the BLR fS mith et al.ll2005[) . This equatorial 
scattering region (ESR) is thus obscured by the torus 
in S2's and has the correct geometry to account for the 
observed alignment of 9y with the radio position angle 
(RPA) in Si's. 

However, it is also evident that Si's as a class exhibit a 
much wider range of opt ical p olarization properties than 
S2's (jSmith et al.ll2062l . i20(M ). In addition to the equa- 
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TABLE 1 

Sample of Polar-scattered Seyfert 1 Galaxies 



Target 


Type 






Ha/H/3 


Pv 




RPA 


Refs. 






(km s^-*-) 


(pc) 




(%) 


n 


n 




Fairall 51 


SB(rs)b Syl 


4114 


55 


5.1 ± 0.6 


3.9 - 7.0 


140 




1 


NGC 4593 


SB(rs)b Syl 


2830 


38 


3.5 ± 0.2 


0.1 - 0.6 


145 




1 


NGC 3227 


SAB(s) pec Syl. 5 


1352 


18 


4.3 ±0.5 


0.6 - 1.6 


125 


170 


1,2 


Mrk 766 


SB(s)a: Syl. 5 


4104 


55 


6.1 ± 0.6 


2.3 - 3.9 


90 


27 


3,4 


Mrk 1239 


E-SO Syl. 5 


5941 


79 


4.7 ±0.8 


4.3 - 6.4 


129 




3 


ESO 323-G077 


SB(1)0"0 Syl. 2 


4417 


59 


5.9 ± 1.7 


2.7 - 4.7 


84 




5 



Note. — Properties of the optically selected sample of polar-scattered Si's. Types have been 
taken from NEEQ. Virgo-infall corrected recessional velocities (vr) have also been taken from NED 
and used to calculate the spatial scale of the _ff5T+NICMOS diffraction limit (0) at 2.0fj,ui (0''2) 
assuming a Hubble constant of 73km s~^ Mpc. We have estimated Ha/H/3 flux ratios from our 
optical spectr opolarimetry of the broad lines (and in the case of ESO 323-G077, from the spectrum 
published by ISchmid et al.l 120031 ) . The optical polarimetry covers the wavelength range 4500 - 
7100A; the range in percentage polarization (red-to-blue) and the average polarization position 
angle over this range are listed. The uncertainties on By are all < 1°. RPA is the radio position 
angle. The mo st up-to-date referen ces for the optical polarimetry (listed first) and radio data (listed 
second) are [ 11 ISmith et al.l i20M). [2] [Mundell et al.i H1995 ). [3] Robinson et al. 2010 (in prep) [4] 
INagar et al.1 lfl999h & [5] ISchmid et al.l I^OM ). 

^http:/ /nedwww. ipac.caltech.edu/ 



torially scattered objects, about 20% exhibit null polar- 
ization, while ^ 25% show characteristics of S2-like po- 
lar scattering. This diversity in the optical polarization 
properties of Seyfert nuclei can be understood if both po- 
lar and equatorial scattering regions are present in all 
Seyferts. The form of the observed polarization is then 
determined largely by the inclination of the torus axis 
to the line-of-sight ([Smith et al.ll2004[ ). See Figure □ for 
a schematic representation of this model. As inclination 
increases from pole-on to edge-on we first see null polar- 
ization Si's, then Si's dominated by equatorial scattering 
(where 9^ is aligned with the radio source axis). With a 
further increase in inclination we see polar-scattered Si's 
(with 0y perpendicular to the radio axis), and finally S2's, 
like NGC 1068, that exhibit polarized broad-lines due to 
polar-scattering. 

In this model, therefore, the polar-scattered Si's oc- 
cupy an orientation regime intermediate between unob- 
scured Si's and a totally obscured S2's. The direct line- 
of-sight to the nucleus passes through the torus at high 
latitudes, and is subject to a modest amount of extinction 
(Ay ~ 1 — 4 magnitudes). The resulting attenuation of 
light from the ESR (within the torus) allows the orthog- 
onally polarized flux from the extended polar scattering 
region (PSR; outside the torus) to dominate the net po- 
larization. In the con t ext of recent clumpy torus models 
(iNenkoya et al.l 120021: iHonig et alj 120061 : iNenkova et alj 
|2008a"b'), an S2 or a polar scattered SI will be observed 
depending on whether or not the line of sight is inter- 
cepted by one or more dense clouds, with the probabil- 
ity of this occurring being a decreasing function of in- 
clination. Unlike in S2's, however, the relatively modest 
line-of-sight extinction in polar scattered Si's does not 
suppress the direct emission from the BLR. 

This simple model can in principle be tested by obtain- 
ing polarimetry observations at longer wavelengths. For 
example, at 2.0/im extinction along direct lines-of-sight 
to the nucleus will be a factor 10 less than at V, sug- 
gesting that it may be possible to detect the polarization 
signature of the ESR at NIR wavelengths, even though 
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Fig. 1.— Refined AGN model of lSmith etHI |(200D showing the 
two-component scattering geometry proposed to explain the optical 
polarization characteristics of Seyfert nuclei. The ESR is modeled 
as a fiared disk that closely surrounds the accretion disk. The PSR 
is modeled as a truncated cone aligned with the axis of the circum- 
nuclear torus. The two regions produce orthogonal polarization. 
We assume that the symmetry axes of the emission disk, both 
the scattering regions and the torus are co-aligned and define the 
principal axis of the system. The observed optical polarization 
depends on orientation: if the line-of-sight lies in region 1 we would 
see a null-polarization SI, in region 2 an equatorially scattered SI, 
in region 3 a polar scattered SI, and in region 4 an S2. 

polar scattering dominates at optical wavelengths. Since 
the scattering geometry dictates that the polarization 
produced by the ESR must be orthogonal to that pro- 
duced by the PSR, the presence of equatorial scattering 
is easily confirmed by comparing the polarization posi- 
tion angle at 2.0/im (02/im) with that measured in the 
optical {9y). Furthermore, the compact ESR will only 
contribute to the polarization of the unresolved nucleus, 
whereas the more extended PSR would be expected to 
dominated any off-nucleus polarization. 

Here we present NICMOS 2.0/im imaging polarimetry 
observations of six polar-scattered Si's, obtained for this 



NIR Polarimetry of "Polar Scattered" Seyfert I's 



3 



TABLE 2 

Details of the Observations 



Target 


Orientation {°) 


S/N 


P 


Fairall 51 


33.4 


281 


1.0 


NGC 4593 


67.0 


193 


1.4 


NGC 3227 


75.6 


106 


1.9 


Mrk 766 


82.6 


57 


2.7 


Mrk 1239 


53.2 


422 


1.0 


ESO 323-G077 


77.0 


525 


0.8 



Note. — Frame orientations for the sample 
(given by the NICMOS "ORIENTAT" header 
parameter) and the theoretically detectable po- 
larization (p) for the measured S/N. The S/N 
was determined from an aperture radius of 7.6 
pixels (0'.'58), corresponding to the unresolved 
nucleus. The values of p are the minimum val- 
ues for which a 3cr detection is possible. 




+/-0.1 
+/-0.2 
+/-0.3 
+/-0.4 
+/-0.5 



400 
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600 



Fig. 2. — Theoretically achievable polarization measurements as 
a function of S/N. The curves correspond to increasing uncertain- 
ties in p (given in the inset). Derived from^Sparks & Axon (1999). 



purpose. The high spatial resolution of HST/NICMOS 
(0'.'2) is desirable as it minimizes cancellation between 
the two orthogonal polarization states and may perhaps 
even resolve the two scattering regions. 

The paper is organized as follows. In §[5] we outline the 
sample selection and describe our observations and data 
reduction. Our approach to the polarimetric analysis of 
the data is described in § |3] and the results are presented 
in § m We discuss these results in the context of the 
two component scattering model in § [S] and present our 
conclusions in § [Bl 

2. OBSERVATIONS AND DATA REDUCTION 

The six Si's selec ted for this s tudy were drawn from 
the list identified bv lSmith et al.l ()2004[ ) as showing clear 
spectropolarimetric signatures of polar scattering. In 
general, polar scattered Si's are characterized by a sys- 
tematic increase in the degree of polarization over the 
optical spectrum, while the polarization position angle 
is approximately constant, without the large excursions 
over the broad Balmer lines that are commonly seen in 
equatorially scattered Si's. For the objects studied here, 
the range in py over the wavelength range ~4500 - 7100A 
and the value of 9^, averaged over the same range, are 
listed for each object in Table [1] along with several other 
properties of interest. 

The observations were designed to mitigate or facili- 
tate removal of several known NICMOS camera 2 (NIC2) 
defects such as residual bad pixels, latent image persis- 
tence, and the inter-pixel response functions (IPRFs). A 
three point spiral dither pattern, with a point spacing 
of S'.'O (40 pixels), was executed a total of four times 
through each polarizer using an exposure time of 10 sec- 
onds (a total of 12 pointings per polarizer). The patterns 
were offset from each other by S'.'O. These observations 
sample the nucleus of each target and minimize the pos- 
sibility of persistence, thoroughly sample the IPRF and 
allow the detection of residual bad pixels. To sample the 
background, 10 second images were taken through each 
polarizer after a 60'.'0 offset. To check for latent image 
persistence, we performed aperture photometry on indi- 
vidual reads in the multiaccum images. Curves of growth 
(count rate vs. time) were constructed and found to be 
linear over the entire exposure. In addition, the residuals 
between image sections both before and after illumina- 



tion by the bright point sources were checked, and found 
to be consistent with the generic noise characteristics of 
NIC2. Therefore, we find no evidence of the latent im- 
ages expected to result from persistence. 

All data were first passed through the calnica pipeline 
with all calibration switches turned on. Accurate polari- 
metric analysis requires that each individual pointing be 
treated separately. Therefore these data were not pro- 
cessed through the calnicb pipeline. The sky images 
were used to subtract the background and residual dark 
current. As the dithered exposures have not been pro- 
cessed as an association (calnicb), they still suffered 
from residual bad pixels and "grot", i.e., pixels with re- 
duced throughput due to particulate contamination on 
the detector. In these cases we used the data quality 
image extensions, and the default "mask- file" to create a 
bad pixel mask that was a combination of hot and cold 
pixels, and pixels affected by grot. These pixels were then 
repaired using a 2D interpolation from the surrounding 
unaffected pixels. 

3. POLARIMETRIC ANALYSIS 

The use of NICMOS as an imaging polarimeter 
has been discussed previously by several authors (e.g., 
iHines et al.l [2OOOI : iBatcheldor et al.l [2001 [200I . Our 

polarization measurements were mad e using the pro- 
cedure s and calibrations described in IBatcheldor et al.l 
(|2009r ). Radial profiles in the Stokes parameters Q and 
U were constructed for each individual pointing using 
the digiphot package within IRAJ0. The average per- 
pixel radial profiles were then computed using an itera- 
tive 2a clipping procedure. In general, rejected profiles 
were affected by uncorrectable bad pixels and "grot", 
and typically amount to 10% of the data for each source. 
However, in the case of ESO 323-G077 three of the 12 
profiles had to be removed. The la dispersion around the 
average for the pointings is used to define the statistical 

uncertainty at each ra dial point. 

Using the recipe of iSparks &: Axoiil ()1999[ ). the final 
Stokes parameter profiles were combined to determine 

^ IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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TABLE 3 

The 2.0/im Polarization Properties of the Sample 



Target 


P2/im(%) 






P2/im(%) 




S9C) 




(0('58) 


(0'.'58) 


{0('58) 


(I'.'S) 


{1"5) 


(1-5) 



Fairall 51 2.2 ± 0.4 139 ± 8 1±8 1.9 ±0.4 148 ± 8 8±8 

NGC 4593 0.4 ±0.6 •■■ ■■■ 1.5 ±0.4 77 ± 6 68 ± 6 

NGC 3227 0.5 ±0.5 ■■■ ■•■ 0.7 ±0.4 102 ± 20 23 ± 20 

Mrk 766 0.6 ±0.7 •■■ ■•■ 2.7 ±0.6 43 ± 5 47 ± 5 

Mrk 1239 0.7 ±0.4 83 ± 9 46 ± 9 0.9 ± 0.3 63 ± 8 66 ± 8 

ESQ 323-G077 1.0 ±0.4 95 ± 10 11 ± 10 0.9 ± 0.3 74 ± 8 10 ± 8 



Note. — Values of P2fim and 62/im for the sample. Values of f?2/xm and S9 
are taken from apertures of radius 0'.'58 and 1'.'5 and are given in degrees in the 
celestial (on sky) reference frame. The values of 59 give the differences between the 
NIR and the average optical polarization position angle (Table [1} • No significant 
polarization was measured within the nuclear apertures of NGC 4593, NGC 3227 
and Mrk 766. 



TABLE 4 
Annular Polarimetry Results 



Target 



Pa{%) 



Fairall 51 
NGC 4593 
NGC 3227 
Mrk 766 
Mrk 1239 
ESQ 323-G077 



0.2 ±0.4 
0.5 ±0.3 
0.4 ±0.1 
1.6 ±0.4 
1.5 ±0.2 
0.9 ±0.2 



41 ± 9 ■ • ■ 36 ± 15 

166 ± 7 ■ ■ ■ 44 ± 27 

24 ± 9 ■ • ■ 19 ± 14 

144 ±6 61 ± 15 81 ± 14 

83 ±8 12 ±18 9 ±16 



Note. — Percentage polarization (pa) and polarization 
position angle 9a at 2fim for a circum-nuclear annulus ex- 
tending from 0'.'58 to 1'.'5. The last two columns give the 
difference in 92p,m between the annular value {9a) and the 
aperture values for radii of 0'.'58 and 1''5, respectively. 



the percentage polarization, position angle and associ- 
ated uncertainties as functions of aperture radius. The 
frame orientations of the observations, which must be 
subtracted from the calculated polarization position an- 
gle in order to determine in the celestial reference 
frame, are given in Table [5] (column 1). We define the 
unresolved nucleus to be the area enclosed by the 2"'^ 
Airy maximum of the HST+NICMOS PSF, which cor- 
responds to 7.6 pixels or 0'.'58. The second column of 
Table[2]gives the integrated S/ N within this nucl ear aper- 
ture for each object. Following Spar ks fc Axonl . we have 
computed the detectable polarization and its associated 
uncertainty as a function of S/N. The polarization de- 
tectable to a given uncertainties is shown for several val- 
ues of the latter in Figure [2j The limiting polarization 
detectable at the 3cr level for the S /N within the nuclear 
apertures of our objects is listed in the third column of 
Table El 

The radial profiles in p2/imand 02/im for each member 
of the sample are presented in Figures [3] and 01 On these 
diagrams, the solid bars indicate the statistical uncer- 
tainty. We note that the statistical errors are correlated, 
as P2/imand 6'2/im are determined from radial integrals 
of the Stokes fluxes. The p2/im panels also indicate the 
upper limit of 0.6% to the inst rumental polariz a tion o f 
NICMOS 2, as determined by Batc heldor et all ()2009[ ). 
The average optical polarization position angle (Table [TJ 



is shown in the ^2/im panels. The NICMOS polarime- 
try is also subject to a systematic error, associated with 
the uncertainties in the parallel transmission coefficients 
tk, which will be the same for all objects. While this can 
shift the values of both p2/im and 02/.tm, the radial trends 
will not be affected. The dashed lines in Figures [3] and |4] 
indicate the allowed ranges of these shifts. Table |3] lists 
the values P2^m and 6'2/im corresponding to aperture of 
radii of 0'.'58 and I'.'S. 

The radial profiles contain information on the spatial 
variation of polarization. However, any extended source 
of polarization will be contaminated by the signal from 
the unresolved nucleus. Therefore, we have also used the 
techniques outlined above to determine the polarization 
within an annular aperture extending to 1'.'5, but exclud- 
ing the inner 0'.'58 of each target. The resulting values of 
P2/im and 6'2/im are given in Table S] 

4. RESULTS 

We have detected significant polarization in all six 
objects, although the radial profiles exhibit diverse 
behavior. At first sight, the radial profiles derived 
from aperture polarimetry suggest that in three objects 
(NGC 3227, Mrk 1239 and ESO 323-G077), p2^in is con- 
sistent with the upper bound on the instrumental po- 
larization. However, the measured values of 02/im differ 
significantly from object to object in the instrumental 
plane, indicating that the measured polarization can- 
not be solely instrumental in origin in all three ob- 
jects. Furthermore, Mrk 1239 exhibits polarization sig- 
nificantly higher than the instrumental upper limit in 
the off-nuclear annulus, as will be discussed in more de- 
tail below. Two of our sample, Fairall 51 and Mrk 1239, 
are also found in NIR polarization study of lBrindle et al.l 
(i 1 990,') . For Fairall 5 1 . our results are consistent with 
iBrindle et al.l 's J & H photo-polarimetry, which was ob- 
tained from a 6'.'0 aperture. The case of Mrk 1239 is 
discussed below. 

Only one source, Fairall 51, exhibits radial polariza- 
tion profiles consistent with a point source. It's p2/im 
and ^2/j,m profil es are similar to those presented by 
iBatcheldor et al.l ()2009[ ) for polarized standard stars, 
which typically have large fluctuations at small radii, due 
to low S/N, but stabilize at constant values of p2/im and 
^2/xm beyond the third Airy ring. In addition, no signifi- 
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Fig. 3. — Radial polarization profiles at 2.0fj,m, derived from aperture polarimetry, for Fairall 51, NGC 4593 and NGC 3227. At 
each point, the polarization is determined from the enclosed Stokes fluxes. The error bars show the statistical uncertainties arising from 
photometric repeatability. The dashed lines show the possible range in values corresponding to the systematic uncertainty on the parallel 
transmission coefficients. [Top Row:] Radial percentage polarizations at 2.0/im. The solid thick line at p2/im = 0.6% is the upper limit to 
the instrumental polarization. [Bottom Row:] 2.0/xm position angles. The thick solid lines mark the optical polarization position angles 
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Fig. 4.— As Fig|3]but for Mrk 766, Mrk 1239 and ESQ 323-G077. 
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cant polarization was detected in the off-nuclear annulus 
for this object (Table HJ. The value of derived 
from aperture polarimetry is the same, to within la, as 
the average optical polarization position angle (Tables [T] 
andEl). 

All 5 of the remaining objects show evidence for ex- 
tended sources of polarization. In fact, the unresolved 
nucleus is not significantly polarized in 3 of these ob- 
jects - NGC 4593, NGC 3227 and Mrk 766. In NGC 4593 
and Mrk 766, the p2/im profiles show a strong increase 
with aperture radius once this exceeds 0'.'58, resulting in 
significant polarizations of, respectively, 1.5 ± 0.4% and 
2.7 ± 0.6% within the 1'.'5 aperture (Table [S]). Both ob- 
jects also exhibit polarization in the off-nuclear annulus 
(ahhough only at the « 2a level in NGC 4593). How- 
ever, in both cases the values of 02^m, whether mea- 
sured from the 1'.'5 aperture or the annulus, differ widely 
from the optical polarization position angle. NGC 3227 
yields only a marginal (« 2a) detection of p2/im in the 
1'.'5 aperture. However, this object is significantly, albeit 
weakly, polarized in the off- nuclear annulus. The values 
of ^2//m measured from the aperture and annulus bracket 
the average optical value (Tables [3] and EJ . 

ESQ 323-G077 is polarized at the level of w 1% in both 
the unresolved nucleus (a 2.5a detection) and in the off- 
nuclear annulus. The values of 02/im measured from the 
two regions are consistent with each other and also with 
the average optical position angle. 

Finally, Mrk 1239 is unique among the objects studied 
here in that it exhibits strong evidence for a change in 
polarization state between the nucleus and the surround- 
ing annulus. The unresolved nucleus is weakly polarized; 
our data yield only sl~ 2a detection. However, the radial 
profile in p2/im shows evidence for a systematic increase 
to « 1% beyond 1'.'4. 6'2/im also varies systematically 
with radius, decreasing from w 80° in the unresolved 
nucleus to « 60° for radii > 1'.'5. The orientation of the 
polarization vector at 2/im therefore differs by w 50 — 70° 
from the average optical value. 

In the 0'.'58 - 1'.'5 annulus, however, the polarization 
is significantly higher, at 1.5 ± 0.2%, than in the nucleus 
and also quite closely aligned (within 15°) with the av- 
erage optical polarization (Tables [1] and S]). Mrk 1239 
is the only object in the sample, for which the annular 
polarization has a position angle significantly different to 
those measured in either the 0'.'58 or 1'.'5 apertures (Ta- 
ble |4]). Perhaps more significant is the comparison with 
ESO323-G077, which is the only other object in which 
we have detected polarization in both the unresolved nu- 
cleus and the annulus. In this object, the 2/im polariza- 
tion in both regions has the same position angle as the 
optical spectrum. In Mrk 1239, on the other hand, there 
is a > 60° change in position angle between the nucleus 
and the circum-nuclear annulus, with the latter having 
approximately the same position angle as measured in 
the optical. 

The systematic change in 02/im with increasing aper- 
ture radius (Fig. |4]) can be understood as the result of 
mixing of the Stokes fiuxes from the unresolved nucleus 
and the surrounding extended source of polarization im- 
plied by the annular polarimetry. As the aperture radius 
increases, and encompasses the PSF, the contribution to 
the integrated Stokes fluxes from the extended source in- 



creases relative to that from the nucleus, causing a grad- 
ual change in 62 um to angles intermediate between the 
intrinsic values for the two sources. In fact, in Fig. |4l 
02^m changes towards the direction of the supplemen- 
tary angle of da- 

Our annular polarization measurements for Mrk 1239 
are consisten t with the J & H photo-polarimetry of 
iBrindle et al.l ()1990[) . This suggests that the extended 
source, not the nucleus, dominates the polarization mea- 
sured within the 6'.'0 aperture used by these authors. 

5. DISCUSSION 
5.1. Predictions of the unified scattering model 

The polar scattered Si's studied here are remarkable 
in that while their total flux optical spectra are typical of 
the general population of Si's, their polarization spectra 
show characteristics similar to those of "hidden broad- 
line" S2's. As previously outlined, these properties can 
be understood if polar-scattered Si's occupy an orienta- 
tion regime (Fig. [1} in which the direct line-of-sight to 
the nucleus follows a relatively transparent path through 
the torus, and is thus subject to modest visual extinction 
{Av 1 - 4 mag). 

In our NICMOS images, the compact ESR should con- 
tribute to the polarization only within the unresolved 
nucleus, whereas the PSR should dominate off-nucleus, 
if it is resolved. If the ESR is the dominant source of 
polarized flux at 2/im, the nuclear point source will have 
a polarization position angle perpendicular to that ob- 
tained from the optical spectropolarimetry. Conversely, 
we would expect any circum-nuclear polarized fiux to 
come from the PSR and hence have a polarization vector 
aligned with that of the optical polarization. 

5.2. Has the equatorial scattering region been detected? 

Our 6 targets exhibit a range of 2/im polarization prop- 
erties, indicating that the origin of the NIR polariza- 
tion is more complex than envisaged in the basic two- 
component scattering model. 

5.2.1. Mrk 1239: equatorial and polar scattering 

This object presents the most compelling case for the 
detection of an ESR. The simplest interpretation of the 
annular polarimetry is that the ESR dominates the polar- 
ization of the unresolved nucleus, whereas at radii >0'.'6, 
where the infiuence of the nuclear PSF is negligible, the 
polarization is dominated by the PSR. The clear detec- 
tion of the ESR signature within 0'.'6 implies that the size 
of this region is less than ^ 200 pc, whereas the annular 
aperture within which the PSR dominates the measured 
polarization corresponds to a linear scale ^ 200 — 600 pc. 
However, it should be noted that the inner radius of the 
annulus was chosen to exclude the nuclear PSF; there is 
no reason to believe it represents a physical boundary 
between the two regions. Indeed, we expect the ESR to 
be located within the torus, with a size limited by the 
dust sublimation radius (~ 0.05 pc for Mrk 1239; see Sec- 
tion I5.3.3|) . It is also likely that the PSR extends much 
further inwards than 200 pc, but contributes a smaller 
fraction of the combined stokes fluxes within the nuclear 
aperture than the ESR. 

Mrk 1239 is a Narrow Line Seyfert 1 galaxy (jGoodrichI 
[1989.) that has a number of remarkable properties. Its 
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optical polarization is among the highest measured in a 
SI, and spectropolarimetry shows that the broad Balmer 
lines are redshifted in polarized flux. This indicates 
that the polar-scattering medium is under going outflow 
m a ^ 1000km s-i wind (' Goodricbl [19891 Robinson et 
al., 2010a, in preparation). It also has an exception- 
ally steep soft X-ray spectru m, which exhibits warm 
absorbers (|Grupe et al.l I2004D. and a strong NIR ex- 
cess peaking at ~ 2^m (jRodrfguez-Ardila &: Mazzalavl 
l200fi[) . A comparison of spectral energy distributions 
leads iRod rfgucz-Ardila & Mazzalay to infer that the 
AGN in Mrkl239 is subject to an intrinsic reddening 
E{B — V) = 0.54, corresponding to a visu al extinction 
Ay ^ 1.7 mag. This is consistent with the iSmith et al.l 
(|2004f) estimate of the extinction along the direct path to 
the nucleus in polar-scattered Si's, and therefore it is en- 
tirely plausible that we are "seeing through" to the ESR 
in the NIR. The NIR excess can be attributed to black 
body emission from hot dust at a tem perature ~ 1200 K 
(|Rodnguez-Ardila fc Mazzalavl[2006l ). 

5.2.2. NGC4593 and Mrk766 

Two other objects, NGC4593 and Mrk766, also ex- 
hibit large differences (59 > 60°) between the 9y and 
6'2/im PA's. However, in both cases the polarization of 
the unresolved nucleus (within '^0'.'6) is consistent with 
zero, the integrated polarization only becoming signifi- 
cant at larger radii. The measured values of P2fim and 
^2/im must therefore be associated with an extended 
scattering region. In these objects, therefore, large val- 
ues of 69 do not imply the detection of the compact 
ESR. In fact, we would expect off- nuclear polarization 
to be due to polar scattering and hence parallel to the 
optical polarization. What then does account for the 
different polarization position angles in the optical and 
NIR? The measured polarization position angle in either 
wave-band could be infiuenced by other sources of polar- 
ization. An obvious source is dichroism due to aligned 
dust grains either in the host galaxy, or our own. This is 
likely to be particularly significant in NGC4593, which 
has a r elatively low op t ical p olarization (0.5%). For ex- 
ample, iHoffman et al.l ()2005[ ) estimate that the optical 
polarization angle could be uncertain by up to 90° due 
to Galactic interstellar polarization (ISP) alone. Thus, 
we have to regard the comparison between the 2/Ltm and 
optical polarization PA's as highly uncertain in this case. 

On the other hand, ISP contamination cannot simi- 
larly explain the observed large 69 observed in Mrk766, 
which has mu ch higher optical polarization (2-4%). 
IHoffman et al.l estimate that Galactic ISP would change 
the optical polarization PA by < 5°. In this case, there- 
fore, we must conclude that polarized flux arising from a 
different scattering geometry makes a substantial contri- 
bution to the off-nucleus NIR polarization. The nature 
of this scattering geometry is unclear, except that the 
average scattering plane must be ~ 45° different from 
that of the region responsible for producing the optical 
polarization. 

5.2.3. Polar scattering 

In Fairall 51 and ESO323-G077 the 2^m polarization 
is closely aligned with the optical position angle, indi- 
cating that polar scattering dominates in the NIR. The 



PSR appears to be unresolved in Fairall 51, implying a 
linear size scale < 150 pc. A simple luminosity scaling of 
the PSR size inferred for Mrkl239 (r oc L^/^) predicts 
that this region should be resolved in Fairall 51; the fact 
that it is not suggests a more compact distribution of 
scatterers in this object. Significant circuni-nuclcar po- 
larization is detected in ESO 323-G077, indicating that 
the polar scattering region is resolved. Evidence for ex- 
tended polarization is also detected in NGC 3227, and 
we similarly attribute this to polar scattering, although 
the value of p2/im is relatively low and 02/xm is not as 
clearly aligned with 9^. The low annular polarization 
is probably the result of geometrical cancellation of the 
polarization vectors. 

5.3. The complex 2^m polarization properties of 
polar- scattered Seyfert 1 's 

Despite having similar optical polarization spectra, our 
small sample of 6 objects exhibits 3 different kinds of be- 
havior at 2/im. In this section, we discuss some possible 
causes of the complexity evident in our results. 

5.3.1. Spatial resolution 

Within the unresolved nucleus it is possible that both 
the equatorial and polar scattering regions contribute to 
the polarized flux at 2^m. The Stokes fluxes from these 
two sources will tend to cancel; if for example, the polar- 
izations are precisely orthogonal, the result will be a net 
polarization either parallel or perpendicular to the sys- 
tem axis, depending on which region produces the higher 
stokes flux. The net polarized flux will be lower than 
would be the case if only one of the scattering regions 
contributes. In principle, this effect could account for the 
weak polarization in the unresolved nuclei of Mrk766, 
NGC 4593 and NGC 3227. Similarly, in Fairall 51 and 
ESO 323-G077, the fact that the PSR appears to dom- 
inate the 2/zm polarization does not necessarily imply 
that the ESR is not present in these objects, simply that 
it produces relatively smaller Stokes fluxes at 2/im. 

If spatial resolution is a key factor in determining the 
observed polarization, we might expect the 2/im polar- 
ization properties to be related to the linear scale (listed 
in Table [T]). However, there is no evidence that this is 
the case, suggesting that the spatial resolution of the 
observations (and consequent mixing of orthogonal po- 
larization states) is not the sole cause of the range of 
observed behavior. Indeed Mrkl239, in which the two 
scattering regions are resolved, is the most distant object 
in our sample. In reality, the relative contributions of the 
ESR and PSR to the net polarization in the unresolved 
nucleus will be determined by multiple variables includ- 
ing the geometry and density of the scattering medium, 
the physical size of the PSR and, as we discuss below, 
the line of sight extinction. 

5.3.2. Extinction 

If the amount of dust extinction along the direct line- 
of-sight to the nucleus is the dominant effect in deter- 
mining whether or not the polarization signature of the 
ESR is actually observed in a given object, we might ex- 
pect a correlation between nuclear reddening and the ob- 
served 2/im polarization properties. Unfortunately, de- 
termining the nuclear reddening in AGN is problematic 
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(|Grandil Il983[ ). The most accessible nebular diagnos- 
tic, the Balmer decrement, does not have a well deter- 
mined intrinsic value for the physical conditions pertain- 
ing in the broad line region. For exampl e, cloud ensemble 
photoionization models computed by iKorista fc GoadI 
|2b04), indicate Ha/H/3 ratios in the range 3.7-4.9. On 
the other hand, recent empirical studies of the broad- 
line Balmer decrement in quasars indicate a mean value 
only slightly great er than the Case B recombination value 
(|Dong et al.l 120081 ). Within our sample, there is a fairly 
wide range in the broad-line Ha/H/3 ratio (3.5 to 6.1, 
Tabled]), but there is no correlation with 59. Moreover, 
the uncertainties are such that there is no significant dif- 
ference in Ha/H/3 between Mrkl239, the object which 
(based on our results) we would expect to be least af- 
fected by reddening, and either Fairall 51 or ES0 323- 
G077, the objects we would expect to be most affected. 
If extinction differences are indeed responsible for the 
observed range in polarization behavior, then these vari- 
ations are masked by the uncertainties in the measure- 
ment of the Balmer decrements, or by variations in the 
intrinsic broad-line Balmer decrements. 

5.3.3. Scattering geometry 

We must also consider the possibility that the effective 
scattering geometry for 2/^m radiation differs from that 
seen by optical photons. In the two component scatter- 
ing model, the ESR is located within the torus, and has 
a flattened, annular structure co-axial with the accretion 
disk. In order to produce the characteristic features of 
the Ha polarization, notably 9{X) variations across the 
line, it must also closely surrou nd the Ha emitting region 
of the disk (jSmith et al.ll2005D . The NIR emission from 
AGN is believed to be dominated by thermal radiation 
from dust heated by the optical- UV continuum of the ac- 
cretion disk; the 2/im continuum coming predominantly 
from the hottest dust near the inner edge of the torus. It 
is unclear if this dust is part of the torus itself , or whether 
it is a distinct component. IMor et all (|2009f ) find that a 
hot black body component (T ~ 1400 X), in addition to 
a clumpy torus, is required to fit the IR spectral energy 
distributions of low redshift quasars, with the black body 
component dominating at wavelengths < 4/im. Similarly, 
hot dust components have been found in several Seyfert 
galaxies, including Mrk 1239 ( as already noted) , Mrk 766 
and NGC4593 (see Table 1 in IRiffel et afl 120091 ). 

Whether or not it is part of the torus, the emitting 
region must be relatively compact; the radius at which 
dust grains have an equilibrium temperature T ^ 1500 K 
is rd ~ 0.4L45^ pc, where L45 is the optical-UV luminos- 
ity in units of 1 0"^^ ergs~^ (jNenkova et al.|[2008al : see also 
IBarvainiill987D . If we take values of XL ° , measured 

. . 5100A' 

from our spectropolarimetric data, to be characteristic 
of the optical-UV luminosity of our low-to-moderate lu- 
minosity Seyferts, then L45 ^3x10^'^ — 2xl0~^ and 
hence Trf 30 — 60 Id. 

The BLR lies within t he torus but also s cales in size 
approx imately as L°-^ (IKaspi et al.l 120051 : IBentz et al.l 
I2009aj) . Using iBentz et aLl 's calibration of the BLR ra- 
dius - optical luminosity scaling relation, we find rBLR ~ 
10 Id for our objects. 

These calculations suggest, therefore, that the 2/im 
dust-emitting region is typically a factor ~ 5 larger than 




Fig. 5. — (a) Polar scattering of NIR emissions from the inner 
edge of the torus. The solid line marks the light path. The solid 
line with double arrows shows the position angle of the polarization 
vector. In this case it is perpendicular to the principal axis due to 
polar scattering. Here, the ESR is a factor ~ 5 smaller than the 
NIR emitting region, (b) Equatorially scattered 2fim emission from 
the outer edge of the accretion disk. In this case, the position angle 
of the polarization vector is aligned with the principal axis, hence 
equatorial scattering is detected. 

the BLR, and hence also the ESR, if the latter closely 
surrounds the BLFiH- The ESR may therefore be much 
smaller than the hot dust NIR emission source and will 
intercept and scatter only a small fraction of the emitted 
flux. On the other hand, the PSR must extend beyond 
the torus and must be at least comparable in size with 
and probably much larger than the NIR-emitting inner 
region. Therefore, polarization of 2/im emission due to 
hot dust located in or near the inner region of the torus 
is likely to be dominated by polar scattering (Fig. [5^). 

However, the outer regions of the accretion disk 
are also expected to emit NIR radiation. Recently, 
iKishimoto et al.l ()2008D reported that in several quasars, 
the NIR continuum in polarized light exhibits an cx 
v^^^ spectrum, consistent with the characteristic long 
wavelength spectrum of the accretion disk. It appears 
that in these objects, it is the accretion disk emission, 
not the torus dust emission that is polarized. As the 
quas a.rs they studied d o not exhibit broad line polariza- 
tion, IKishimoto et al.l argue that the scattering region 
is located well outside the NIR emitting region of the 
accretion disk but within the radius of the BLR. A sim- 
ple black body accretion disk model indicates that the 
disk will be cool enough (^ 1500 K) to emit at 2/im at a 
radius r2/im - 7 x lO^Rg{L/ LEddV^HM,/lO'' Mq)-^/^ , 
where LEdd is the Eddington luminosity, Rg — 2GM,/c^, 
and we have assumed a radiative efficiency of 10%. 
For the lower luminosity Seyferts studied here, black 
hole masses estimat ed either by reverberation map- 
ping (|Peterson et al.l [2004: Bentz et al. 2009B), or by 
the virial method (e.g., iVestcrgaard & Peterson 120061 : 
iGreene fc Hoi 120051 see also Peterson et al. 2004) lie in 
the range M, ~ 0.2 - 4 x 10 ''M0. Therefore, given that 
L/LEdd < 1, the simple disk model suggests that the 
2/im emission comes from radii < 10 Id, a region compa- 
rable in size with Balmer line emitting BLR. In general, 
therefore, the NIR emitting region of the disk lies close 
to but within the ESR and accordingly we expect that 
2/im accretion disk emission should be mainly subject to 
equatorial scattering, producing polarization parallel to 
the disk axis (and hence perpendicular to the observed 

Note, however, that in a NIR rever beration ma pping study 
of 4 nearby Seyferts, including NGC 3227, ISuganuma et al. (200® 
find that the NIR emission region is only a factor ~ 2 larger than 
the BLR. 
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optical polarization, Fig. [SJd)- 

These arguments suggest that the observed 2/im polar- 
ization may in general depend on the relative importance 
of scattering of torus emission by the PSR, or scattering 
of accretion disk emission by the ESR. As our objects are, 
by selection, dominated by polar scattering in the opti- 
cal, and as dust emission dominates the AGN continuum 
at wavelengths > 2/im, it might be expected that po- 
lar scattering will dominate the 2/im polarization even 
if the the direct line of sight to the ESR is essentially 
transparent. This appears to be the case in Fairall 51 
and ESO 323-G077. On the other hand, the two scatter- 
ing regions are apparently resolved in Mrkl239. We can 
plausibly attribute the polarization of the unresolved nu- 
cleus to equatorial scattering of the accretion disk emis- 
sion, whereas the polarization detected ofF-nucleus is due 
to polar scattering of the hot dust emission. Thus, while 
the latter component might dominate the total 2/im flux 
of the unresolved nucleus, it does not contribute to the 
polarization in the nucleus. 

The inner walls of the torus itself may also scatter 
2^m emission from the disk and/or hot dust components. 
Scattering of light from a central point source by the 
inner walls of an optically thick torus can produce po- 
larization either parallel to or perpendicular to the axis, 
depending on the torus opening angle and inclination 
(IKartid 119951 : ICorbett et"al][2000l: iGoosmann fc GaskeUl 
120071 ) As the dust sublimation radius is unresolved, this 
polarization would only contribute to the net polariza- 
tion of the unresolved nucleus. However, such a contri- 
bution is unlikely to be significant. If the NIR emission 
is dominated by hot dust in or near the torus opening, 
i.e., approximately co-spatial with the scattering region, 
there will be a wide range in scattering angles leading 
to cancellation of the Stokes fluxes and a low net polar- 
iz ation. On the other hand, th e monte-carlo simulations 
of IGoosmann fc Gaskelll (|2007f) indicate that the polar- 
ized flux produced by torus scattering of disk emission, 
i.e., approximately a central point source, will be much 
smaller than the polarized fluxes produced by the PSR or 
ESR, for any plausible electron scattering optical depth. 

5.4. Wavelength dependence of polarization degree 

A key property of the polar scattered Seyferts is that 
the optical continuum polarization increases strongly to 
the blue. This trend continues into the NIR, in the 
sense that the 2/im polarization of the unresolved nu- 
cleus is significantly lower, in most objects, than that 
measured at ^ 0.5/(m. The exceptions are NGC4593 
and NGG 3227, both of which have relatively low optical 
polarization and comparably low polarization at 2/im. 
The wavelength dependence of p{X) can be attributed, 
at optical wavelengths, to dilution by the reddened, di- 
rect AGN (accretion disk) continuum, with a contribu- 
tion from unpolarized stellar emission (which is required 
to account for the local increase s in p over the broad 
emission lines; iSmith et al.ll200^ . An old stellar popu- 
lation would also dilute the polarization at 2//m. How- 
ever, it seems likely that the dominant effect in the un- 
resolved nucleus is dilution by hot dust emission. As 
already noted in Section I5.3.3[ black body components 
attributed to hot dust emission have been identified in 
several Seyfert galaxies, including Mrk 1239, Mrk 766 and 
NGG 4593. In the latter two cases, in particular, our 



aperture polarimetry reveals behavior consistent with the 
presence of a diluting component associated with the nu- 
cleus: the 2/tm polarization increases significantly with 
radius, while the nucleus itself is weakly polarized at a 
level < 0.6%, the upper limit for instrumental polariza- 
tion. 

We note in addition that the wavelength dependence 
of p(A) also depends on the composition of the scattering 
medium. Dust scattering may well dominate in the PSR 
for which, if in the Rayleigh regime (grain sizes < A), 
the scattering cross-section oc A"**. Unpolarized diluting 
sources (e.g., starlight; accretion disk and hot dust emis- 
sion) will therefore have a much greater effect in the NIR 
than in the optical band. 

6. GONCLUSIONS 

We have measured the nuclear polarization at 2.0/im in 
six "polar scattered" SI galaxies in an effort to confirm 
the presence of a compact equatorial scattering region 
in these objects. This in turn would support the two- 
component scattering model for Seyfert galaxies. We 
find evidence for both equatorial and polar scattering 
in Mrk 1239, which exhibits the expected signatures of 
equatorial scattering in the unresolved nucleus, and po- 
lar scattering off-nucleus on linear scales of > 100 pc. 
In this source, therefore, the compact ESR is revealed 
in NIR scattered light and is clearly resolved from the 
PSR, which extends beyond the torus. A further impli- 
cation is that the source of the scattered 2/im emission 
in the unresolved n ucleus is the accretion d isk, as in the 
quasars studied bv iKishimoto et all ()2008[) , rather than 
torus hot dust emission. 

The remaining objects exhibit a variety of behavior, 
suggesting that in general the origin of the 2/im polar- 
ization is more complex than envisaged in the simple 
two component scattering model. In Fairall 51 and 
ESO 323-G077, the nuclear polarization is evidently 
dominated by polar scattering, even at 2/im. This 
cannot be taken as evidence of the absence of an ESR, 
it simply implies that the Stokes fluxes produced by 
this region are smaller than those produced by polar 
scattering. This could be the result of higher extinction 
to the nucleus, a smaller covering factor, or it could 
be that the scattered 2/tm radiation in these objects 
comes predominantly from hot dust in or near the inner 
part of the torus which 'see's the PSR rather than 
the ESR. In both NGG 4593 and Mrk 766, aperture 
polarimetry shows an increase in the 2/im polarization 
with radial distance from the nucleus (where, indeed, it 
is consistent with zero) . Gonceivably, this behavior could 
be the result of cancellation between the orthogonal 
polarization components of the two scattering regions 
within the unresolved nucleus. However, in both cases, 
there is a large difference between the 2/tm and optical 
polarization position angles off nucleus, where polar 
scattering is expected to dominate. This may not be 
significant for NGG 4593, in which contamination by 
interstellar polarization is likely to be important, but 
is perplexing in the case of Mrk 766, which is strongly 
polarized in the optical. The nucleus of NGG 3227 is 
also not significantly polarized at 2/tm. However, in 
this object the surrounding PSR is resolved in deeper 
imaging polarimetry, to be presented elsewhere. 
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